AD-A108  757  COASTAL  ENGINEERING  RESEARCH  CENTER  FORT  BEL VOIR  VA  F/G  8/10 

PRODUCTS  FROM  TWO  COMPUTER  PROGRAMS  WHICH  PROCESS  DI8ITAL  BATHT— ETC(U) 
OCT  81  BE  HERCHENRODER 

UNCLASSIFIED  CERC-CETA-81-13 _ NL _ 


* 


MkrxuOOPY  RESOLUTION  TE ST  CHART 


mr ADA108757 


CETA  81-13 


Products  From  Two  Computer  Programs 
Which  Process  Digital  Bathymetric  Data 


by 

Barry  E.  Herchenroder 

COASTAL  ENGINEERING  TECHNICAL  AID  NO.  81-13 

OCTOBER  1981 


A 


Approved  for  public  release; 
distribution  unlimited. 


U.S.  ARMY,  CORPS  OF  ENGINEERS 

COASTAL  ENGINEERING 
RESEARCH  CENTER 

Kingman  Building 
Fort  Belvoir,  Va.  22060 


Reprint  or  republication  of  any  of  this  material 
shall  give  appropriate  credit  to  the  U.S*  Army  Coastal 
Engineering  Research  Center. 

Limited  free  distribution  within  the  United  States 
of  single  copies  of  this  publication  has  been  made  by 
this  Center.  Additional  copies  are  available  from: 

National  Technical  Information  Service 
ATTN:  Operations  Division 
5285  Port  Roual  Road 
Springfield,  Virginia  22161 

The  findings  in  this  report  are  not  to  be  construed 
as  an  official  Department  of  the  Army  position  unless  so 
designated  by  other  authorized  documents. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  r*h*n  Dele  Entered) 


REPORT  DOCUMENTATION  PAGE 


.  REPORT  NUMBER 

CETA  81-13 


ft.  TITLE  (and  Subtltla) 

PRODUCTS  FROM  TWO  COMPUTER  PROGRAMS  WHICH 
PROCESS  DIGITAL  BATHYMETRIC  DATA 


7.  AUTHORfaJ 

Barry  E.  Herchenroder 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Department  of  the  Army 

Coastal  Enginering  Research  Center  (CERRE-CO) 
Kingman  Building,  Fort  Belvoir,  Virginia  22060 


II.  CONTROLLING  OFFICE  NAME  AND  AODRESS 

Department  of  the  Army 

Coastal  Engineering  Research  Center 

Kingman  Building,  Fort  Belvoir,  Virginia  22060 


MONITORING  AGENCY  NAME  ft  ADDRESSfff  different  from  Controlling  Office) 


1«.  DISTRIBUTION  STATEMENT  (of  thla  Report) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT'S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  ft  PERIOD  COVEREO 

Coastal  Engineering 
Technical  Aid 


6.  PERFORMING  ORG.  REPORT  NUMBER 


ft.  CONTRACT  OR  GRANT  NUMBERfa) 


10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


A31592 


12.  report  date 


October  1981 


IS.  NUMBER  OF  PAGES 


IS.  SECURITY  CLASS,  (of  thlo  report) 


UNCLASSIFIED 


ISa.  OECLASSIFICATION/DOWNGRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  ebatrect  entered  in  Block  20,  it  different  from  Report) 


19  KFY  WORDS  (Continue  on  rove  tee  elde  It  nec  eeeery  mid  Identity  by  block  number) 


Computer  programs 


Digital  bathymetric  data 


20  AB£T*RACT  (Continue  on  reweree  elde  II  neceeeery  end  Identify  by  block  mmaber) 

A  description  is  given  of  products  from  two  computer  programs  which  process 
digital  bathymetric  data.  One  program  generates  regularly  spaced  bathymetric 
data  from  irregularly  spaced  data.  The  other  uses  regularly  spaced  data  to 
determine  and  draw  contours.  A  large  set  of  irregularly  spaced  bathymetric 
data  available  on  magnetic  tape  for  U.S.  coastal  regions  is  also  described. 
Fxamples  of  output  from  each  program  are  displayed  for  two  coastal  areas. 


DO 


EDITION  OF  t  NOV  ftft  IS  OBSOI  FTE 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dete  Entered) 


PREFACE 


This  report  describes  products  from  two  computer  programs ,  and  a  set  of 
digital  bathymetric  data  for  U.S.  coastal  regions.  The  programs  and  data  are 
useful  for  generating  regularly  spaced  bathymetric  data  and  determining  and 
drawing  contours.  The  work  was  carried  out  under  the  waves  and  coastal  flood¬ 
ing  research  program  of  the  U.S.  Army  Coastal  Engineering  Research  Center 
(CERC) . 

The  present  report  was  prepared  by  Dr.  Barry  E.  Herchenroder,  under  the 
general  supervision  of  Dr.  C.L.  Vincent,  Chief,  Coastal  Oceanography  Branch, 
Research  Division. 

Comments  on  this  publication  are  invited* 


Approved  for  publication  in  accordance  with  Public  Law  166,  79th  Congress, 
approved  31  July  1945,  as  supplemented  by  Public  Law  172,  88th  Congress, 
approved  7  November  1963. 


3 


CONTENTS 

CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI) 

I  INTRODUCTION  . 

II  NOS  DIGITAL  DATA  . 

Ill  PROGRAM  ZGRID1 . 

IV  PROGRAM  CONTOR . 

V  LIMITATIONS . 

1.  ZGRID1  . 

2.  CONTOR  . 

VI  SUMMARY . 

LITERATURE  CITED  . 


Page 

6 

7 

7 

8 

8 

8 

15 

15 

15 

16 


FIGURES 


1  Index  (1°  square  areas)  to  hydrographic  data  (except  those  of  the 

Alaskan  coast)  ,  1930  to  1965  only .  9 

2  Index  (1°  square  areas)  to  hydrographic  data  off  the  Alaskan  coast, 

1930  to  1965  only . 10 

3  ZGRID1- produced  soundings  offshore  of  Great  Egg  Harbor  Inlet, 

New  Jersey,  in  a  0,03°  square  area . 11 


4  ZGRID1- produced  bed  elevations  near  CERC’s  FRF  at  Duck,  North  Carolina  -  12 


5  ZGRID1- produced  regularly  spaced  depth  values  for  the  same  area  off 

of  Great  Egg  Harbor  Inlet,  New  Jersey,  covered  by  Figure  3 . 13 

6  CONTOR- produced  depth  contours  for  the  same  area  off  of  Great  Harbor 

Inlet,  New  Jersey,  covered  by  Figures  3  and  5 . 14 

7  CONTOR- produced  depth  contours  for  the  same  region  near  CERC's  FRF 

covered  by  Figure  4 . 14 


5 


•  i 


i  i  i  i  Trfaaik 


CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


U.S.  customary  units 
metric  (SI)  units  as 

of  measurement 
follows: 

used  in 

this  report  can  be  converted  to 

Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x 

o 

i 

u> 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

lTo  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  formula:  C  -  (5/9)  (F  -32). 

To  obtain  Kelvin  (K)  readings,  use  formula:  K  -  (5/9)  (F  -32)  +  273.15. 


PRODUCTS  FROM  TWO  COMPUTER  PROGRAMS  WHICH 
PROCESS  DIGITAL  BATHYMETRIC  DATA 


by 

Barry  E .  Herahenroder 
I.  INTRODUCTION 

Regularly  spaced  bathymetric  data  are  needed  by  the  U.S.  Army,  Corps  of 
Engineers  for  input  to  various  numerical  models  which  describe  coastal  phenom¬ 
ena  such  as  storm  surge  and  wave  refraction.  Bathymetric  contours  produced 
from  such  data  are  also  sometimes  needed  to  help  locate  important  bathymetric 
features  or  changes  in  a  project  area.  To  generate  regularly  spaced  bathy¬ 
metric  data,  an  adequate  set  of  raw  irregularly  spaced  data  must  be  available 
and  a  method  to  interpolate-extrapolate  this  data  to  regularly  spaced  locations 
is  required.  This  report  describes  a  large  set  of  bathymetric  data  available 
on  magnetic  tape  and  products  from  a  computer  program  which  generates  regularly 
spaced  data  from  the  raw  data;  products  from  a  program  which  generates  contours 
from  the  regularly  spaced  data  are  also  described. 

One  source  of  raw  bathymetric  data  in  an  area  affected  by  a  coastal  project 
may  be  survey  data  taken  as  part  of  the  project.  If  there  are  no  such  data  or 
if  the  project  survey  data  must  be  supplemented,  a  major  additional  data  source 
is  the  National  Ocean  Survey  (NOS)  digital  hydrographic  data  set  available  from 
the  National  Geophysical  and  Solar-Terrestrial  Data  Center  (NGSDC)  of  the 
National  Oceanic  and  Atmospheric  Administration  (NOAA) .  Data  from  both  data 
sources  normally  have  an  irregular  spacing. 

Interpolating-extrapolating  irregularly  spaced  data  to  regularly  spaced 
locations  is  onerous  when  done  by  hand.  Drawing  depth  or  bed  elevation  con¬ 
tours  by  hand  is  also  time  consuming  and  tedious.  To  simplify  these  tasks  at 
the  Coastal  Engineering  Research  Center  (CERC) ,  two  computer  programs,  ZGRIDl 
(program  No.  720X6RICG0)  and  CONTOR  (program  No.  720X6R19X0),  have  been  devel¬ 
oped  or  modified.  ZGRIDl  interpolates  or  extrapolates  and  CONTOR  determines 
and  draws  contours.  The  two  programs  are  available  with  documentation  from  the 
Engineering  Computer  Programs  Library  at  U.S.  Army  Engineer  Waterways  Experiment 
Station,  Vicksburg,  Mississippi.  Both  ZGRIDl  and  CONTOR  require  large,  fast 
computers,  e.g.,  CDC  6600,  CDC  CYBER176,  CRAY-1,  IBM  3033,  and  use  standard 
CALCOMP  plotting  commands. 

II.  NOS  DIGITAL  DATA 

A  few  years  ago,  NOS  began  digitizing  U.S.  coastal  hydrographic  data  from 
about  3,200  survey  sheets.  Most  of  the  resulting  digital  data  have  been 
archived  with  the  NGSDC;  they  are  available  on  magnetic  tape  for  a  fee  by  writ¬ 
ing  or  phoning  National  Geophysical  and  Solar-Terrestrial  Data  Center,  Code 
D621,  EDS/NOAA,  Boulder,  Colorado  80302,  phone  (303)  499-1000,  extension  6338 
(FTS  373-6338).  A  detailed  description  of  the  available  digital  data,  as  well 
as  pricing  information,  is  given  in  National  Geophysical  and  Solar-Terrestrial 
Data  Center  (1976,  1979a,  1979b,  1980)  and  in  Lawrence  (1977).  The  digital 
data  are  irregularly  spaced  in  latitude  and  longitude  and  consist  of  about  97 
percent  depth  soundings  and  3  percent  bottom  characteristics  (e.g.,  clay,  mud, 
soft,  etc.)  and  dangers  to  navigation  (e.g.,  wrecks,  pilings,  etc.)  The  1° 
square  areas  along  the  continental  United  States  and  Alaskan  coasts  for  which 
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digital  data  were  available  as  of  January  1981  are  indicated  in  Figures  1  and 
2.  Data  from  the  squares  in  Figures  1  and  2  were  digitized  from  surveys  con¬ 
ducted  between  1930  and  1965. 

III.  PROGRAM  ZGRID1 

ZGRID1  is  a  modification  of  the  program  supplied  by  NGSDC  to  perform  two- 
dimensional  interpolation-extrapolation.  Irregularly  spaced  bathymetry  data 
from  the  NOS  digital  tapes  or  from  independently  obtained  surveys  are  read  by 
ZGRID1.  ZGRIDl  then  determines  the  coordinates  of  each  sounding  on  a  map  pro¬ 
jection  preselected  by  the  user  and  interpolates-extrapolates  this  data  to 
regularly  spaced  locations  on  the  projection. 

As  an  option  of  ZGRIDl,  a  plot  of  the  irregularly  spaced  soundings  can  be 
made.  Figures  3  and  4  are  examples  of  such  plots.  Figure  3  shows  the  sound¬ 
ings  offshore  of  Great  Egg  Harbor  Inlet,  New  Jersey,  on  a  0.03°  (1.8  minute) 
latitude-longitude  square.  Soundings  were  obtained  from  the  NGSDC  tape  covering 
the  1°  latitude-longitude  square  39°  to  40°  N.  latitude,  74°  to  75°  W.  longitude. 
Figure  4  consists  of  bed  elevations  near  CERC's  Field  Research  Facility  (FRF)  at 
Duck,  North  Carolina.  The  FRF  bed  elevations  were  derived  from  a  survey  made  by 
CERC  in  September  1978. 

The  estimated  bed  elevation  or  sounding  value  at  each  regularly  spaced  loca¬ 
tion  can  be  written  on  a  magnetic  disk  or  tape  for  subsequent  use.  At  the 
user’s  discretion,  the  regularly  spaced  values  can  be  printed  on  a  lineprinter 
or  plotted.  Figure  5  is  a  plot  of  the  estimated  regularly  spaced  depth  values 
for  the  same  region  covered  by  Figure  3.  For  this  example,  there  are  30  regu¬ 
larly  spaced  (on  the  map  projection)  points  in  both  the  x  and  y  directions  for 
a  total  of  900  regularly  spaced  points  in  the  whole  region. 

IV.  PROGRAM  CONTOR 

C0NT0R  determines  and  draws  contours  of  any  variable  from  a  two-dimensional 
array  of  variable  values.  These  values  are  assumed  to  be  located  at  regularly 
spaced  grid  points.  CONTOR  was  developed  at  the  Carnegie-Mellon  University, 
Pittsburgh,  Pennsylvania,  by  Robinson  and  Scarton  (1972).  When  used  to  contour 
bed  elevation  or  water  depth  fields,  CONTOR  is  usually  run  in  conjunction  with 
ZGRIDl.  ZGRIDl  produces  the  depth  or  bed  elevation  values  at  regularly  spaced 
locations  needed  by  CONTOR.  The  contours  to  be  plotted  are  specified  by  the 
user . 

Figures  6  and  7  are  examples  of  contour  plots  done  by  CONTOR.  Figure  6 
corresponds  to  the  region  off  Great  Egg  Harbor  Inlet  covered  by  Figures  3  and 
5.  Figure  7  corresponds  to  the  FRF  region  covered  by  Figure  4.  Figure  6  was 
produced  from  the  ZGRIDl  computed  values  in  Figure  5.  Figure  7  was  produced 
from  ZGRIDl  computed  regularly  spaced  values  (not  shown)  located  82.02  feet 
(25  meters)  apart  in  both  x  and  y  directions. 

V.  LIMITATIONS 

Limited  experience  indicates  that  the  data  on  the  NGSDC  tapes  contain  a 
few  gross  errors  in  most  0.25°  by  0.25°  squares.  One  satisfactory  way  to 
detect  these  errors  is  through  the  use  of  contour  plots.  Gross  data  errors 
an:  easily  located  on  these  plots  since  they  appear  as  small  "bull’s-eyes.1' 
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Index  (1°  square  areas)  to  hydrographic  data  off  the 
Alaskan  coast,  1930  to  1965  only,  available  on  tape  from 
NGSDC  (National  Geophysical  and  Solar-Terrestrial  Data 
Center,  1980). 
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Figure  3.  ZGRID1- produced  soundings  offshore  of  Great  Egg  Harbor  Inlet, 
New  Jersey,  in  a  0.03°  square  area.  Axis  annotations  are  N. 
latitudes  (positive  numbers)  and  W.  longitudes  (negative 
numbers) ,  Depths  are  in  tenths  of  meters  referenced  to  MLW. 
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Figure  4.  ZGRID1- produced  bed  elevations  near  CERC's  FRF  at  Duck,  North  Carolina, 
The  CERC  pier  has  been  drawn  in  by  hand.  Axis  annotations  are  in  meters 
Bed  elevations  are  in  tenths  of  meters  reference  to  the  1929  National 
Geodetic  Vertical  Datum.  Negative  elevations  are  above  the  datum. 


Figure 


1.  ZGRID1 . 


ZGRID1  sometimes  gives  poor  results  near  a  land-water  boundary  or  an  outer 
computational  boundary  where  there  are  too  few  irregularly  spaced  raw  data 
points.  Where  there  is  a  significant  topographic  trend,  ZGRID1  can  also  give 
misleading  results  if  the  number  of  raw  data  points  is  too  small. 

2.  CONTOR. 

When  a  bathymetric  feature  has  a  size  comparable  to  the  grid  spacing,  con¬ 
tours  which  go  near  or  through  the  feature  may  cross  or  have  unrealistic  loops 
or  wiggles.  To  minimize  this  contour  "noise,”  the  regularly  spaced  depth  or 
bed  elevations  produced  by  ZGRTT>1  may  have  to  be  smoothed. 

VI .  SUMMARY 

This  report  gives  general  descriptions  of  products  from  two  computer  pro¬ 
grams  (ZGRID1  and  CONTOR)  which  process  digital  bathymetric  data.  ZGRIDl  gen¬ 
erates  regularly  spaced  bathymetric  data  from  Irregularly  spaced  data,  and 
CONTOR  uses  regularly  spaced  data  to  determine  and  draw  contours.  A  large  set 
of  bathymetric  data  available  on  magnetic  tape  for  U.S.  coastal  regions  is 
also  described.  Examples  of  program  output  are  presented  for  two  areas,  one 
near  Great  Egg  Harbor  Inlet  and  the  other  near  CERC’s  FRF. 
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